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The total r hadronic width can be accurately calculated using analyt- 
icity and the operator product expansion. The result turns out to be very 
sensitive to the value of Q!s(m^), providing a precise determination of the 
strong coupling constant. The theoretical description of this observable is 
updated, including the recently computed O^a"^) contributions. The exper- 
imental determination of asim^) and its actual uncertainties are discussed. 

PACS numbers: 12.38.-t, 12.38.Bx, 12.38.Qk 

1. Introduction 

The inclusive hadronic decay width of the r lepton provides one of the 
most precise measurements of the strong coupling [l:-_7j. Moreover, the 
comparison of as(m^) with Os determinations at higher energies constitutes 
the most accurate test of asymptotic freedom, successfully confirming the 
predicted running of the QCD coupling at the four-loop level. 

The calculation of the 0(af) contribution [8] has triggered a renewed 
theoretical interest on the as(m^) determination, since it allows to push 
the accuracy to the four-loop level [8lll3j. However, the recent theoretical 
analyses slightly disagree on the final result, giving rise to a range of different 
values for as{m'^). The differences among these results, shown in Table [H 
are too large compared with the claimed accuracy and originate in 

the different inputs or theoretical procedures which have been adopted. 

In the following I try to clarify the reasons behind these numerical dis- 
crepancies and reassess the actual uncertainties of the r decay determination 
of as- Using all present experimental and theoretical knowledge, I derive 
the value 

a^(m2) = 0.342 ±0.012. (1) 



* Presented at the Flavianet topical workshop on Low energy constraints on extensions 
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Reference 


Method 


Sp 




Baikov et al. [8j 
Davier et al. [9] 
Beneke-Jamin |10| 
Maltman-Yavin [11] 

Menke [12] 
Caprini-Fischer |13j 


CIPT, FOPT 

CIPT 
BSR + FOPT 
PWM + CIPT 
CIPT, FOPT 
BSR + CIPT 


0.1998 ±0.0043 
0.2066 ±0.0070 
0.2042 ± 0.0050 

0.2042 ± 0.0050 
0.2042 ± 0.0050 


0.332 ±0.016 
0.344 ± 0.009 
0.316 ±0.006 

0.321 ±0.013 

n Q42 + 0-011 

U.OZU _ Q QQg 



Table 1. C'(af) determinations of asim^). 



2. Theoretical framework 

The inclusive character of the total r hadronic width renders possible 
an accurate calculation of the ratio [IHS] 

p - r[T- ^ hadrons (7)] p , p , p /9^ 

T[t- Vre~Ve(n)\ 

The theoretical analysis involves the two-point correlation functions for the 
left-handed quark currents L^ - = ■0j7'*(l — 75)'0i (i, j = u, d, s): 

ni;^(g) ^ i J d'xe^'^^ (o|r(Lf^.(x)Lr,(o)t)|o) 

= (-5'"^^' ± g'^Q") + gMg- . (3) 

Using analyticity, Rr can be written as a contour integral in the complex 
s-plane running counter-clockwise around the circle |s| = m^: 



= 4) f 1 ^ 



i±24 n(o+i)(5)-24nW(s) 



(4) 

where n('^)(s) = |V^dpnirf^(s) ± l^wsP nii^(s). This expression requires 
the correlators only for complex s of order , which is significantly larger 
than the scale associated with non-perturbative effects. Using the Opera- 
tor Product Expansion (OPE) to evaluate the contour integral, Rr can be 
expressed as an expansion in powers of 1/m^. The uncertainties associated 
with the use of the OPE near the time-like axis are heavily suppressed by 
the presence in dH) of a double zero at s = [HEJE]. 

The theoretical prediction for the Cabibbo-allowed decay width can be 
written as [1] 

Rry+A = Nc \Vud? Sew {1 ± 5p ± 5np} , (5) 
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where Nc = 3 is the number of quark colours and Sew = 1-0201 ± 0.0003 
contains the electroweak radiative corrections jl5H17j . The dominant cor- 
rection (~ 20%) is the perturbative QCD contribution in the massless-quark 
hmit 5p, which is akeady known to 0(af) [T|l8] . Quark mass effects p]ll8pi9j 
are tiny for the Cabibbo-allowed current and amount to a neghgible correc- 
tion smaher than 10~^ pillOj. 

Non-perturbative contributions are suppressed by six powers of the r 
mass fl] and, therefore, are very small. Their numerical size has been deter- 
mined from the invariant-mass distribution of the final hadrons in r decay, 
through the study of weighted integrals |20j which can be calculated theo- 
retically in the same way as R-r- The predicted suppression [IJ of the non- 
perturbative corrections has been confirmed by ALEPH [21j . CLEO [22] 
and OPAL [23j- The most recent analysis [6] gives 

5np = -0.0059 ±0.0014. (6) 

The measured values of the r lifetime and leptonic branching ratios im- 
ply Rr = 3.640 ± 0.010 [9]. Subtracting the Cabibbo-suppressed contribu- 
tion Rr^s = 0.1615 ± 0.0040 [9], one obtains Rr,v+A = 3.479 ± 0.011. Using 
\Vud\ = 0.97418 ± 0.00027 [21] and ([6]), the pure perturbative contribution 
to R-r is determined to be: 

5p = 0.2038 ± 0.0040 . (7) 



3. Perturbative contribution to R,- 

In the chiral limit (m„ = rud = rag = 0), the vector and axial- vector 
currents are conserved. This implies sn*^'^)(s) = 0; therefore, only the cor- 
relator n(''"'"^)(s) contributes to Eq. The result is more conveniently 
expressed in terms of the logarithmic derivative of the two-point correlation 
function of the vector (axial) current, n(s) = ^ n(*'"^^)(s), which satisfies an 
homogeneous renormalization-group equation: 

With the choice of renormalization scale fj? = Q'^ = —s all logarithmic 
corrections, proportional to powers of log {—s/fj,'^), have been summed into 
the running coupling. The Kn coefficients are known to order a^. For 
nj = 3 quark flavours, one has [8l[25H29]: Kq = Ki = 1, K2 = 1.63982, 
= 6.37101 and xf^ = 49.07570. 
The perturbative component of Rr is given by 



6p = J2 KnA^^\as), 

n=l 



(9) 



4 



Loops 




A(2)(a.) 


A(3)(a,) 






Sp 


1 


0.13247 


0.01570 


0.001698 


0.000169 


0.0000154 


0.1773 


2 


0.13523 


0.01575 


0.001629 


0.000151 


0.0000124 


0.1788 


3 


0.13540 


0.01565 


0.001597 


0.000145 


0.0000115 


0.1784 


4 


0.13529 


0.01557 


0.001579 


0.000142 


0.0000112 


0.1779 



Table 2. Exact results for A'-'^'>{ag) (n < 5) at different /3-function approximations, 
and corresponding values of Sp = X]ri=i A''"'\cts), for Q!s(m^)/7r = 0.1. 



where the functions [2] 

A('-)(as) ^ — i - ( ^iiz^V (i_2— + 2—-—] (10) 
27ri J\s\=m2 s \ n J \ ml ml j ^ ' 

are contour integrals in the complex plane which only depend on the strong 
coupling. Using the exact solution (up to unknown /3n>4 contributions) for 
ois{s) given by the renormalization-group /3-function equation, they can be 
numerically computed with a very high accuracy ^2j. 

Table [2] gives the numerical values for A^'^\as) {n < 5) obtained at the 
one-, two-, three- and four-loop approximations (i.e. /3n>i = 0, /3ra>2 = 0, 
/^n>3 = and /3n,>4 = 0, respectively), together with the corresponding re- 
sults for 5p = Y^n=i KnA^'^\as), taking as(m^)/7r = 0.1. The perturbative 
convergence is very good. The error induced by the truncation of the (3 func- 
tion at fourth order can be conservatively estimated through the variation 
of the results at five loops, assuming = ±/3l/f33 = =f443, i.e. a geometric 
growth of the /3 function. Higher-order contributions to the Adler function 
D{Q^) will be taken into account adding the fifth-order term K^A^^\as) 
with = 275ib400. Including, moreover, the 5-loop variation with changes 
of the renormalization scale in the range iJ? G [0.5,1.5], one gets the final 
result 6p = 0.1810 ± 0.0045^^ ± Q.mUp^ ± 0.0013^ = 0.1810 ± 0.0049 for 
Usim^)/!! = 0.1. 

Adopting this very conservative procedure, the experimental value oi 6p 
given in Eq. ([7]) implies the strong coupling determination quoted in Eq. ([T]) . 

4. Discussion on previous 0{q.^) determinations of Q:s(m^) 

The integrals A^'^\as) can be expanded in powers of Or = as(m^)/7r, 
^H(a,) = + 0(a^+i). One recovers in this way the naive perturbative 
expansion (5p = X]n=i (-^n + 5") = ^n=i '''nO-r This approximation 
is known as fixed-order perturbation theory (FOPT), while the improved 
expression ([9]), keeping the non-expanded values for the integrals A^"-\as), 
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is usually called contour-improved perturbation theory (CIPT) [2 ll30j . FOPT 
gives rise to a pathological non-convergent series (its radius of convergence 
is slightly smaller than the physical value of Or), because the long running 
of as{s) along the circle |s| = generates very large gn coefficients, which 
depend on Km<n and ^rn<n ffi = 0, 52 = 3.56, 53 = 19.99, 54 = 
78.00, g5 = 307.78. These corrections are much larger than the original Kn 
contributions and lead to values of as(m^) smaller than ([T]). FOPT suffers 
from a huge renormalization-scale dependence f2]. As shown in a recent 
detailed analysis |12j . the actual FOPT uncertainties are much larger than 
usually estimated. Once this is taken properly into account, the FOPT 
results are consistent with CIPT, but their huge uncertainties make them 
irrelevant. 

An ad-hoc model of higher-order coefficients for the Adler function has 
been recently advocated |1U) . The model mixes three different types of 
renormalons (n = —1, 2 and 3) plus a linear polynomial, trying to enforce a 
cancelation of the Kn and gn coefficients in order to get a better behaviour 
of the FOPT series. It contains 5 free parameters which are determined 
by the known coefficients of the Adler series and the chosen value of K^. 
Making a Borel summation of the full renormalon series (BSR), one gets a 
sizeable positive contribution to 6p from higher orders, implying a smaller 
value for as(m^) independently of the adopted FOPT [TO] or CIPT [13] 
procedure. This model constitutes an interesting example of possible higher- 
order corrections, making apparent that the associated uncertainties have to 
be carefully estimated. However, it cannot be used to determine the physical 
value of as, because the result is model dependent. Different assumptions 
about the unknown Kn coefficients would lead to different central values for 
as- 

The determination of as(m^) in Eq. ([T|) agrees with Refs. [6l[T2] and 
with the CIPT value in Ref. It is somewhat larger than the (CIPT) 
result extracted in Ref. [11] from pinched-weight moments (PWM) of the 
hadronic distribution, which seems to correspond to a very different value of 
6]yp- After evolution up to the scale Mz ^Mj; the strong coupling in Eq. ([T|) 
decreases to 

a,(M|) = 0.1213 ±0.0014, (11) 
in excellent agreement with the direct measurements at the Z peak. 
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